Critical amplitude ratios of the Baxter- Wu model 
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A Monte Carlo simulation study of the critical and off-critical behavior of the Baxter- Wu model, 
which belongs to the universality class of the 4-state Potts model, was performed. We estimate the 
critical temperature window using known analytical results for the specific heat and magnetization. 
This helps us to extract reliable values of universal combinations of critical amplitudes with reason- 
able accuracy. Comparisons with approximate analytical predictions and other numerical results 
are discussed. 



I. INTRODUCTION 



One of the central results of the theory of phase tran- 
sitions and critical phenomena is the formulation of the 
universality hypothesis [l|, y] . According to the theory, 
all systems with the same dimensionality, the same sym- 
metry of the ordered phase and the same number of order 
parameters are described by the same set of critical expo- 
nents at the critical point. Additionally, thermodynamic 
functions vary with temperature in such a way that some 
combinations of their amplitudes take the same values 
for all systems within a universality class 3- For many 
systems, critical exponents are known by approximate 
methods (field theoretical perturbation theory, series ex- 
pansions, Monte Carlo simulations) and have been de- 
rived exactly in some cases, mainly for two-dimensional 
models [^-Q ■ The extensive research in the last 50 years 
gives strong support for the universality hypothesis in 
the context of the critical exponents [5g]. At the same 
time, the issue of critical amplitude ratios was checked 
only for a few models [Ij]. 

Special interest in the properties of universality classes 
derives from cases in which the singular behavior is com- 
plicated by logarithmic corrections [15|. There are some 
systems which belong to the same universality class but 
whose off-critical behavior may be modified by logarith- 
mic corrections. This is the case for the universality 
class of the two-dimensional 4-state Potts model. The 
model, which gives the name to this universality class, 
the 4-state Potts model [16] contains logarithmic correc- 
tions to the critical behavior [17-20] of thermodynamic 
quantities. For example, the free energy F and the mag- 
netization M in the ordered phase behave at criticality 

as \n,M 



Fir) « A4p|rr"[ln(-r)]"' 
Mi-r) « i34p(-r)^[ln(-r)]A, 



(1) 
(2) 



where a=2/3, /?=1/12, a;=-l, ^==-1/8, and 
T—\~Tc/T is a measure of the distance of the tempera- 
ture T to the critical temperature Tc. 

Another two-dimensional model, which belongs to the 
universality class of the 4-state Potts model, is the 
Baxter- Wu model Wa. For this model it is known ex- 



actly [21| - |24| that close to the critical temperature, the 
free energy and the magnetization in the ordered phase 
behave as 



M{~t) « Bbw(-T)^, 



(3) 

(4) 



with the same values for the critical exponents as in ^ 
and ([2]), but without logarithmic corrections to the sin- 
gular behavior. 

A possible explanation for the difference in the off- 
critical behavior of the Baxter- Wu and 4-state Potts 
models is that for some unknown reason the coefficient 
behind the log arithmic correction is zero for the Baxter- 
Wu model [ij, [3 , although the leading critical behavior 
for both models seems to be described by the same fixed 
point in the renormalization-group space [25] . 

Critical amplitudes are not universal and one should 
not expect that the free-energy amplitudes ^bw of the 
Baxter- Wu model and ^4p of the 4-state Potts model 
are equal. At the same time, some combinations of the 
amplitudes are universal. For example, the ratio of the 
free-energy amplitude Abw(+) in the high-temperature 
phase of the Baxter- Wu model to the amplitude ^bw(— ) 
in the low-temperature phase equals unity, as well as the 
ratio of the corresponding amplitudes of the 4-state Potts 
model ^4p(-|-)/yl4p(— )=1. This is a consequence of the 
duality relation for the free energy of the two models [2J] . 

The ratio of the susceptibility amplitudes 
7^x=^4p(+)/^4p(-)=^bw(+)/^bw(-), on the other 
hand, is not known exactly and different approximations 
(analytical, series expansions, Monte Carlo data) from 
different groups are not coherent and differ substantially 
(for a recent discussion, see the papers [20, [2a|). It is 
well known that logarithmic corrections, if they exist, 
complicate the critical behavior and render the analysis 
of the critical behavior extremely difficult if possible at 
all, and the determination of critical exponents from 
Monte Carlo (MC) and series expansions (SE) became 
non-trivial. In fact, the determination of critical expo- 
nents and corrections to scaling may lead to indecisive 
conclusions 27]. 

The purpose of the present paper is to estimate numer- 
ically critical amplitudes of the Baxter- Wu model, which 



is free of logarithmic corrections, and to compare univer- 
sal amplitude ratios with the ones available for the 4-state 
Potts model. It should be emphasized that critical am- 
plitude ratio universality is a non-local property of the 
renormalization group phase space whereas the critical 
exponents characterize its behavior only in the vicinity 
of the corresponding fixed point. 

We use the traditional Metropolis MC algorithm to 
simulate the Baxter- Wu model, and analyze the magne- 
tization and polarization in the ordered phase, and the 
energy, specific heat and magnetic susceptibility in both 
phases. We estimate the accuracy of our data by compar- 
ing the magnetization, polarization, energy and specific 
heat to available exact results. 

In particular, we use the known results for the energy, 
specific heat, magnetization and polarization in order to 
estimate the critical temperature window, i.e., the range 
of temperatures in which the system on a finite lattice 
behaves to a very good approximation as on an infi- 
nite one. This allows us to estimate the critical ratio 
of the susceptibility amplitudes in the high- and low- 
temperature phase with good accuracy, TZ^ — 3.9 ± 0.1. 
The analytical estimate of this ratio, obtained by Delfino 
and Cardy 28] using some approximation of the ex- 
act scattering field theory [20|, is only slightly larger, 
TZ^ — 4.013. Delfino and Grinza [3^ obtained practically 
the same value TZ^ = 4.02 using the same approximation 
for the Ashkin- Teller model with parameters which corre- 
spond to the 4-state Potts model universality class. A re- 
cent analysis of the MC and SE data for the 4-state Potts 
model gives an amplitude ratio in the range of about 
6.5(4) [20i|2ll3i|. At the same time, the values for the 
universal combination of amplitudes R^ = aA{^T/B'^ 
in the low-temperature phase reported in Ref. [2^ and 
Refs. 0,[2i,IIl| are 0.00508 and 0.0052(2), respectively 
{Aq and Bq are the specific-heat and magnetization am- 
plitudes). This is in perfect agreement with our estimate 
0.00517(7) we present here for the Baxter- Wu model. 

In the rest of the paper we present the details of 
our analysis of the amplitude ratios for the Baxter- Wu 
model. In Section lU and Hill we first give an overview of 
known analytical results and previous numerical simula- 
tions which support them. We then discuss in Section ITVl 
details of our simulation algorithm realization, including 
the special choice of the lattice, averaging, etc. Section IVl 
presents the details of our critical amplitude estimation, 
and the discussion in Section IVII summarizes our results 
and touches on some open questions. 



A. Model 

The Baxter- Wu model is defined on a triangular lat- 
tice, with spins ai = ±1 located at the vertices. The 
three spins forming a triangular face are coupled with a 
strength J, and the Hamiltonian reads 



"H = - J ^ (TiO-jCrfe , 



(5) 



faces 



where the summation extends over all triangular faces of 
the lattice, both pointing up and down. 



B. Self-duality and critical temperature 

The model is self-dual as found by Wood and Grif- 
fiths 32 1 and Merlini and Gruber [3j], who applied the 
Kramers- Wannier construction developed for the square- 
lattice Ising model, and showed that both models, Ising 
on the square lattice and Baxter- Wu on the triangular 
lattice, share the same self-dual temperature (see also 
the book 2J]). The partition function reads 



^-E 



exp 



13 ^ (JtCTjak 

faces 



(6) 



where /3 — J /ksT. The dimensionless free energy per 
site is 

f/keT = HP) = - lim N-^ In Z , (7) 

where N is the number of lattice sites. It satisfies the 
duality relation 



where 



V'(/3)=V(/3*)+ln(sinh2;3* 



sinh2;3*sinh2/3 = 1 



(8) 



(9) 



This is precisely the duality relation of the square lattice 
Ising model. The argument of Kramers- Wannier can be 
applied: if there exists just one critical point, then it 
must occur when /3 = Pc — P*, where 



sinh2/3c = l, /?c = ln(V2 + l)/2. 



(10) 



II. MODEL AND EXACT RESULTS 



C. Ground-state symmetry 



In this section we summarize those known analytical 
results for the Baxter- Wu model which will be used in 
the lattice construction, algorithm realization and data 
analysis. 



Let aA, fJB, CTc denote all the spins on the A, B, 
C sub-lattices, respectively. Any total configuration 
(cr^, (TBjCrc) of spins has the same energy as three oth- 
ers. These configurations can be obtained by fiipping all 



spins on two of the sub-lattices. The spin configurations 
can be grouped in equal-energy sets of four: 

(— cr^jCrs,— CTc) , (— CTA, — CBjCTc) • (H) 

The ground state is thus four-fold degenerate: one ferro- 
magnetic state with magnetization (per site) rn = 1 and 
three ferrimagnetic states with m = 1/3. 



D. Exact solution and critical behavior 

Baxter and Wu solved the model at the critical temper- 
ature exactly [2l| - [24| and found that the critical value of 
the energy e = {'H)/N is eg = — y^J J|, the specific heat 
C = de/dT diverges at T^ as [HI, ^ 



Coc|i|-'/': 



(12) 



and the critical behavior of the magnetization for r < 
is M 



moc|t|i/i2 



(13) 



where here the reduced temperature is defined as t = 
(T — Tc)/Tc- The critical exponents thus take the values 
a=2/3 and /3=1/12, which is the two-dimensional 4-state 
Potts model universality class. 



The magnetization reads [36[ (i < 0) 

m = |t| A (Bo + Bi|<|i + B2\t\ + Bj,\t\^ + 0(|t|i)) , 

(21) 



where 



Bo = 25 (ln(%/2 + l))'' =1 



283 264 709. 



(22) 



1 



Bi = ^ (ln(^/2 + l))* =-0.589 829 210... ,(23) 

28 ^ ' 



2 = ^^(ln(V2 + l))'^(4-^/21n(y2 + l)) 
= 0.073 615 269... , (24) 

1 \ — 

— ('ln(%/2 + 1)) '' ^ 0.542 208 469 6 ... . (25) 



B, = 



Some time later, Baxter 24] found an elegant form 
for the free energy in terms of an infinite sum of the 
ratio of some polynomials. The magnetization (and also 
the polarization defined below in Eqs. p4p and ^51) ) can 
be expressed in terms of an infinite product of ratios of 
polynomials. These expressions may be used to obtain 
the energy, specific heat, magnetization and polarization 
with any desired accuracy. The summary of analytical 
results is presented in this section in a form which is 
most suitable for the analysis of the data discussed here. 



E. Joyce's results for C and M 

Joyce established analytic properties of the free energy 
per spin and found the explicit form 35i] 



/ 



2|J| 



InA(M) , 



where 
1 



1 



1 14: 16m(1-m)2 



A(u) (i + uy 



iFi 



2'6'3' (1-l-u)'' 



(14) 



(15) 



with 2F1 denoting the hypergeometric function and u — 
exp(— 4| Jj/fc^T) such that the critical point value in this 
variable is Uc = 3 — 2v2. 

The specific heat critical behavior follows as 



ks 
where 

v4n 



Aoltr 3 + ^1 + A2t\t\-T^ + Aaltl 3 + 0(t) , (16) 



(ln(\/2- 



= 0.187 787 867. 



(17) 



A, = -i(ln(V2 + i: 
A2 = 4fln(\/2 + l) 



27 
= 0.295 775 490. 



A. 



ln(\/2 + i; 



= -0.388 409 700... ,(18) 

('7\/21n(\/2 + l)-4 

(19) 

= 0.396 723182.... (20) 



III. PREVIOUS NUMERICAL RESULTS 

A. Monte Carlo study of the critical behavior of 
the Baxter- Wu model 

In course of their simulation study of the infiuence of 
quenched impurities on the critical behavior, Novotny 
and Landau jST,] performed for comparison also first 
Monte Carlo analyses of the pure Baxter- Wu model. 
They defined the magnetic order parameter m as the 
root-mean-square (rms) of the magnetization on the 
three sub-lattices. The polarization order parameter p 
was defined as the rms average of the two-spin correla- 
tion functions for the nearest-neighbor spins between dif- 
ferent sub-lattices. They pointed out that the rms rather 
than the the absolute value of the sum of the sub-lattice 
magnetizations was used because the susceptibilities cal- 
culated using the rms definition gave results that were in 
closer agreement with high- and low-temperature series- 
expansion results. 

Finite-size scaling plots of the order parameters, the 
magnetization niL^™-/^ and polarization pL''p/'^, as func- 
tions of tL^/^ are consistent with /3m=/3p=/3=l/12 and 
1^=2/3. In the symmetric phase (T > T^) those plots 
show in a log-log representation a decay with slope 
— 7/12, approaching asymptotic values at very large val- 
ues of iL^''^. The energy, specific heat, and susceptibility 
also behave according to the Baxter- Wu predictions. 

They simulated lattices with linear size L up to 
66 and estimated magnetic susceptibility amplitudes. 



rbw(+)=0.03(2) and rbw(-)=0.010(5), and polariza- 
tion amplitudes, rP^(+)=0.06(3) and rP^(-)=0.04(2). 
Therefore, the ratio of amplitudes is about 1.5 — 3 with 
an uncertainty of about 50 per cent. 



B. Monte Carlo renormalization-group study of 
the Baxter- Wu model 



In their Monte Carlo renormalization-group (MCRG) 
analysis, Novotny, Landau, and Swendsen [38| did not 
find any evidence for logarithmic corrections in the 
Baxter- Wu model on lattices with sizes 21 x 21, 42 x 42, 
and 147 x 147. For the RG eigenvalues they estimated 
2/T=l-48 ± 0.03, y/,=1.875 ± 0.003, and results for 2/3 
were consistent with the Barber ansatz |39i] j/3 = 7/8, 
based on the rcnormalization group analysis. The latter 
exponent is the correction-to-scaling magnetic exponent 
2/^2 = 2 — Xa^ (see, Refs. [3, [3| and the discussion in 
Ref. [20^), which produces the correction-to-scaling expo- 
nent 2/3 in the specific heat ([T^ and magnetization (I^TI) . 



C. Conformal invariance studies of the Baxter- Wu 
model and a related site-coloring problem 

The operator content of the Baxter- Wu model has been 
discussed by Alcaraz and Xavier [40, |4l| , who extended 
the original Bethe ansatz solution of the site-coloring 
problem and solved numerically the corresponding equa- 
tions of the transfer matrix for toroidal boundary con- 
ditions. They found that the latter problem has the 
same operator content as the 4-state Potts model. The 
correction-to-scaling effects seem, however, to correspond 
to different perturbations of the fixed point of the rcnor- 
malization group. The authors of Refs. [40, |41[ claim 
that the correction-to-scaling terms contain only integer 
powers, like those in the Ising model. 



IV. SIMULATIONS 
A. Lattice choice 

Genuine critical behavior can be observed only in the 
thermodynamic limit, when the system approaches in- 
finite size. Simulations, on the other hand, are always 
performed on finite lattices Wi\. Finiteness of the lat- 
tice leads to such effects as scaling of thermodynamic 
quantities with the lattice size at criticality [i^, and the 
shift [43] of the pseudo-critical point [5g . 

An additional source for systematic deviations in the 
simulations is due to the approximation of the aspect 
ratio which is often overlooked. The central idea is to 
choose such a form of the finite lattice and boundary 
conditions for which the number of sites and number of 
bonds would be the same in either direction, and on both 




FIG. 1: (Color online) Top row: Three elements generating 
the lattice with sites R, G, and B (circle, box, and triangle) 
and bonds G (solid line), B (dashed line), and R (dotted line), 
attached correspondingly. Bottom: The smallest lattice hav- 
ing 3 sites, 3 bonds, and zero sum of site colors and bond 
colors in any of the three lattice directions under periodic 
boundary conditions, L=3. 



the prime and its dual lattice [59|. This is the property 
of the model on the infinite lattice. 

We construct the triangular lattice by using the three 
elements shown in the top row of Fig. [T] We associate 
one of the three colors with each site and bond, and use 
the convention that the sum of R (red), G (green), and 
B (blue) vanishes, or, equivalently, equals to W (white) 
color. The rule of construction is that a site of the same 
color should never be a neighbor, and the same condition 
should hold for the bonds. This is the natural partition 
of the Baxter- Wu lattice on the three sub-lattices. We 
denote in Fig. [1] a R-site with a circle, a G-site with a 
box, and a B-site with a triangle. G-bonds are denoted 
in Fig. [1] with a solid line, B-bonds with a dashed line, 
and R-bonds with a dotted line. 

The minimal lattice size which follows the above men- 
tioned requirements is shown at the bottom of Fig. [TJ It 
contains exactly three sites of different colors and three 
bonds of different colors in any of the three lattice di- 
rections with periodic boundary conditions. The sum of 
site colors and bond colors is white along any lattice di- 
rection (by definition, R + G + B = W = Q). This 3x3 
lattice can be used as an elementary building block for 
constructing larger lattices. An example of the next size 
of the lattice, the 6x6 lattice, is shown in Fig. [51 

There are the three sub-lattices for the Baxter- Wu 
model labeled R, G, and B, therefore the total num- 
ber of sites (and bonds) should be a multiple of three. 
The bonds between each two sub-lattices form a hexag- 
onal lattice: for removed sites R it is composed by the 




FIG. 2: (Color online) The next smallest 6x6 lattice with 
6 sites, 6 bonds, and "zero" or "white" colors in any lattice 
direction under periodic boundary conditions, L=6. 



C. Metropolis algorithm 

To update the spin configurations with the Metropolis 
algorithm [4g| we calculate the local energy of a spin at 
position (j, fc), 

ej,k — ^o'j,k (cj^fc-iCj+i^fc + crj+i^kCj+i,k+i 

+ o-j+iM+if^jM+i + o-j,k+icrj-i,k (26) 

+ 0-j-l,kO'j-l,k-l + 0-j-l,k-l<7j,k-l) ■ 

If Cj.fc > 0, we flip spin cFj.k- Otherwise we flip the spin 
only if exp(2/3ej^fc) is not less than a uniformly distributed 
random number G (0, 1]. 



D. Averaging over the ensemble 

In simulations, the specific heat can be found from 



hexagons formed by the red bonds icB, for removed sites 
G by the hexagons formed by the green bonds Ljxb , and 
for removed sites B by the hexagons formed by the blue 
bonds -Lrg- 

Such construction keeps the symmetry of the Baxter- 
Wu model. In addition, the Baxter- Wu model is self- 
dual and our choice of the lattice construction keeps self- 
duality not only in the thermodynamic limit, but also 
for any finite size of the lattice. This in turn minimizes 
the possible influence of the approximation for the aspect 
ratio H. 



B. Choice of the algorithm 

It is well known that Monte Carlo simulations of the 
Baxter- Wu model experience strong finite-size effects and 
an application of the traditional Metropolis algorithm be- 
comes costly due to the critical slowing down. Novotny 
and Evertz [43] proposed some time ago a cluster algo- 
rithm for the Baxter- Wu model. Recently, this algorithm 
has been extended to the generalized self-dual Baxter- Wu 
model [4a] . The main idea is that the Hamiltonian is in- 
variant under the transformation where all spins on the 
two sub-lattices are changed. The algorithm fixes the 
spins on one sub-lattice and builds up clusters of corre- 
lated spins on the remaining two sub-lattices (see Fig. [2]) . 
It is, however, not obvious that such clusters which live 
on the subspace of possible configurations will percolate 
at the critical point of the Baxter- Wu model. Indeed, we 
found that the percolation point of such clusters seems to 
be shifted a little bit to lower temperatures. As a result 
this leads to the shift of the curves for some observables. 
The same effect was found recently for the behavior of 
Fortuin-Kasteleyn clusters in the Z4 spin model \A7l . De- 
spite the slowing-down problem, we therefore resorted in 
our simulations to the traditional Metropolis update al- 
gorithm. 



C^N- 



1 



((^ 



(27) 



where the energy per site is calculated as 

T7 2^ 0-j^k<7j + l.k+li<^j+l,k + fj.fc+l) (28) 



N 

and N = L'^ denotes the number of sites. 

Similarly, the reduced magnetic susceptibility in the 
low-temperature phase, x_, can be obtained from 



kBTx- = N {{m') - {mf) , 



(29) 



where the magnetization m is computed as the sum of 
the magnetization over the three sub-lattices. 



E^ 



(30) 



with the magnetization per site, m^, of the sub-lattice i 
given by 



Li 



(31) 



In the analysis, we actually calculated the magnetiza- 
tion in an alternative way as 13 7] 






(32) 



since it leads to more accurate results as was already 
mentioned in Ref. (33] and will be discussed in more detail 
below. 

The magnetic susceptibility in the high-temperature 
phase is computed as 



kBTx+ ^N{m^) 



(33) 




o 




FIG. 3: Specific lieat in the high-temperature phase for sys- 
tem sizes L = 162 (open boxes), L — 243 (open circles), 
L = 363 (stars), and the four-term approximation Eqs. (|16p - 
(|20[1 to the exact solution (solid line). Error bars are of the 
size of the data symbols. 



The polarization p per site follows as the sum of the 
polarizations 



Pi 



J^^'^3A(^],k-i + o-j+i.fe+i + ctj-im) (34) 



between two sub-lattices j and k: 

3 



P = 



^E 



Pi 



(35) 



We compute the average of a quantity A (where A is 
one of e, e^, e^, m, m^, m'^, nis, m'^, m^, p, p^, p"^) as a 
sum over iVav steps. 



(■'^''^Y.-' 



N, 



Before averaging, we equilibrate the system with iVrciax 
MC steps. Typically, A^rciax = 10^ - 10^ and A^av = 
10^ - 10^. 



E. Dual reduced temperatures 

We compute thermodynamic quantities at the reduced 
temperature values t and r* connected via the duality 
relation ^ which can be written in the form 



FIG. 4: Specific heat in the low-temperature phase. Symbols 
and curves are the same as in Fig. [S] 



The reduced temperatures r and r* are defined as 
T-Tr 



T 
T* -Tc 



= 1 - T,/T, r > , (37) 

= l-Tc/T\ T*<0. (38) 



tanh /3* — e 



-2/3 



(36) 



Due to ([55t the reduced temperatures are related by 

r = 1 + -^ ln{tanh[/3,(l - r*)]} , r* < . (39) 

In some formulas we also employ the reduced tem- 
perature t = '^TfT'" as in Joyce's papers |3a, [Sg and in 
Eqs. ^, dH]), which is related tor by iw t+t'^+0{t'^). 

V. RESULTS 

In this section we employ natural units in which J — 
fcs = 1 and first define the temperature region window 
in which we would fit our data. We use the exact knowl- 
edge of the specific-heat behavior for that purpose. We 
demonstrate how reliable the fits to the data are. 



A. Temperature region window 

The specific-heat data in Figs. [3] and |4] exhibit strong 
finite-size effects close to the critical temperature, that is 
when the reduced temperature t approaches zero. This 
is particularly pronounced in the high-temperature phase 
shown in Fig. [31 where one can see visible deviations of 
the MC data sets from the exact solution for very small 
t. For temperatures t > 0.003, however, the shown data 
sets coincide with each other and with the exact result. 
For that values of temperature, the correlation length 
becomes smaller than the system size and the relation 



S,o oc t '^ <^ L holds better for larger reduced tempera- 
tures. 

For large (absolute) values of reduced temperature, the 
computational data become more and more close to the 
exact values. At the same time, the solid line, which 
represents an approximation to the exact solution, starts 
to diverge from the computational data, because higher- 
order correction-to-scaling terms, which are not included 
in the approximation, become more and more important 
for larger reduced temperatures. So, the critical tempera- 
ture window is bounded for smaller reduced temperatures 
\t\ by finite-size effects and for larger reduced tempera- 
tures |f I by the neglected correction-to-scaling terms in 
the analysis. 

We have to stress that the left boundary of the tem- 
perature window is determined clearly by the nature of 
the phase transition - it is the temperature at which two 
length scales coincide: the (randomly) fixed system size 
and the temperature dependent correlation length. The 
right boundary of the temperature window is not fixed 
by any physical reason. It depends on the correction- 
to-scaling variation with temperature and the number of 
terms considered in Eq. ([T51) . Here some conventions are 
necessary. For example, we can define the right bound- 
ary as the temperature up to which the first two terms of 
the full correction set are important. Of course, one can 
also choose one term or three terms. For some systems or 
for some particular value of a tuning parameter, it may 
happen that the first correction-to-scaling term is close 
to zero, so that in some wide temperature region the sys- 
tem would behave as an infinite one, i.e., the correction- 
to-scaling terms are not importaiit (see the discussion of 
such an extended scaling in Ref. |49|). 

In the case of the Baxter- Wu model, correction-to- 
scaling terms are not small. The power of correction 
terms decays slowly with the exponent 1/3 [see Eq. p^ ]. 
The amplitudes of correction terms do not depend on 
any parameter. So, we have to choose some convention. 
Using the known powers and amplitudes of correction- 
to-scaling terms for both the specific heat and magne- 
tization, we can estimate the right edge of the tempera- 
ture window as that (reduced) temperature for which the 
relative deviation of the truncated expansion (J16p . de- 
noted by Cjv(|t|fl) with N terms, from the exact value of 
C{\t\B) is smaller than some e, \CN{\t\R)IC{\t\R) - 1| < 
e. Fixing some value of e, say e ~ 0.001 (deviation less 
than one tens of per cent), we can estimate the value of 
\t\R as a function of the number of correction-to-scaling 
terms N we want to include in the analysis. 

For the data analysis, we will use a combined set of 
data obtained for the system size L — 363 when |i| < 0.02 
and for the size L = 243 for larger reduced temperatures 
|i|, if not mentioned otherwise. Both data sets are com- 
puted with A^reiax = 10^ MC steps for relaxation and 
A^av = 10^ for averaging. 



B. Specific heat 

We define effective amplitudes of the specific heat by 

AM±){t)^C{t)\t\^/\ (40) 

where the argument "±" refers to the high-temperature 
(reduced temperature t is positive) and low-temperature 
phase (reduced temperature t is negative), respectively. 
Figures [5] and [S] show our MC data for the effective 
specific-heat amplitude in the high- and low-temperature 



TABLE I: Results of fits to the MC data for the effective 
specific-heat amplitude in the high-temperature phase. 



fit window 


M+) 


^i(+) 


M+) 


^3(+) 


exact (Ref. [35]) 


0.1877... 


-0.3884 . . 


. 0.2957. . 


0.3967... 


0.001<f<0.76 


0.1878(8) 


-0.42(3) 


0.91(7) 


-0.29(5) 


0.001<f<0.32 


0.1879(12) 
0.1862(8) 


-0.43(5) 
-0.32(2) 


0.93(17) 
0.58(3) 


-0.31(14) 



0.003<f<0.76 


0.1865(8) 
0.1827(6) 


-0.38(3) 
-0.24(1) 


0.82(6) 
0.46(1) 


-0.24(4) 



0.003<i<0.32 


0.1852(12) 
0.1849(6) 


-0.31(6) 
-0.30(1) 


0.59(18) 
0.54(2) 


-0.04(15) 



0.001<f<0.025 


0.1893(22) 
0.1846(11) 


-0.52(16) 
-0.16(2) 


1.2(5) 







0.001<f<0.015 


0.1860(15) 


-0.20(4) 








0.001<f<0.01 


0.1870(23) 


-0.24(7) 









TABLE II: Resuhs of fits to the MC data for the effective 
specific-heat amplitude in the low-temperature phase. 



fit window Aq{—) 


^i(-) 


M-) M-) 


exact (Ref. [35]) 0.1877... 


-0.3884 . . . 


0.2957. . . 0.3967 . . . 


0.001<|ii<0.4 0.1895(6) 
0.1972(11) 


-0.457(28) 
-0.863(22) 


0.615(88) 1.05(7) 
-0.679(30) 



0.003< 



O.OOK 



0.003< 



O.OOK 



0.003< 



O.OOK 



0.003< 



O.OOK 



<0.4 0.1899(7) -0.45(29) 0.565(88) 1.01(7) 



<0.2 



0.1879(7) 
0.1932(2) 



-0.351(40) 
-0.718(24) 



1.03(15) 1.47(15) 
0.411(41) 



<0.2 0.1874(7) -0.321(39) 1.31(15) 1.55(14) 
0.1945(9) -0.749(24) -0.461(40) 



<0.1 



0.1903(6) 
0.1889(3) 



-0.587(22) 
-0.529(4) 



0.122(45) 




<0.1 



0.1908(7) 
0.1888(3) 



-0.605(23) 
-0.528(4) 



0.156(46) 




<0.025 



0.1888(10) 
0.1895(4) 



-0.482(72) 
-0.541(9) 



0.19(23) 




<0.025 



0.1875(13) 
0.1897(4) 



-0.398(82) 
-0.544(8) 



0.44(25) 




<0.01 0.1885(11) -0.511(30) 








0.20 




FIG. 5: Effective amplitude of the specific heat in the high- 
temperature phase. Symbols show combined MC data for L — 
243 (for t > 0.02) and 363 (for t < 0.02), and the sohd hue 
is the four-term approximation Eqs. H16p - (|20[) to the exact 
solution. 



FIG. 6: Effective amplitude of the specific heat in the low- 
temperature phase. Symbols and curves are the same as in 
Fig. El 



phase, respectively, together with the four-term approxi- 
mation Eqs. p^ -(pn )) to the exact solution, which is seen 
to be in a fairly good agreement with the numerical data. 
We fit the following expression for the effective ampli- 
tude of the specific heat to the data: 

^eff (±)(t) - ^o(±) + Ai{±)\t\i + A2i±)t + A3i±)\t\^ 

(41) 
by varying the temperature region window. Results of 
the fits are presented in Tables HI and HIl The data clearly 
support in both phases the specific-heat amplitude value 
Aq — 0.1877... rather well. The first correction-to- 
scaling amplitude Ai = —0.3884 is estimated less well, 
but is still precise enough to give support for the theo- 
retically expected value 2/3 of the correction-to-scaling 
exponent. Higher correction-to-scaling terms could not 
be estimated from the available MC data set. 



C. Energy 




FIG. 7: Effective amplitude of the energy in the low- 
temperature phase. Symbols are MC data and solid line is 
the four-term series expansion (|43p. (|44|l . 



The specific-heat amplitudes are connected with the 
energy amplitudes. Let us define effective energy ampli- 
tudes as 



eeff(+)(t) = (e+(t)-eo)t-^, 
eoff(-)(i) = {e-{t)~eo)\t\--^ 



(42) 



In the vicinity of the critical point, they can be expanded 
as dl 

eeff(±)(i) = Ei+ E2\t\i + Est + Ei\t\^ + 0(|t|i) (43) 
with coefficients 

El = 1.278 376 401... , 



E2 = -0.881371587. 
E3 = 0.503 377 046.. 
E4 = 0.540143 046.. 



(44) 



Figure [7] shows the effective amplitude of the energy in 
the low-temperature phase together with the four-term 
series expansion (|43)) . p4)) . Fits to the energy amplitudes 
are shown in Table Hill which clearly support the first two 
coefficients in (H^ . Fits to the effective energy amplitude 
in the high-temperature phase look very similar. 



D. Magnetization 

From here on we follow the usual convention in the 
magnetic sector and use the reduced temperature r as 
independent variable in figures and fits. 

The magnetization may be estimated using two meth- 
ods, defined by Eqs. ([30]) and ^. Figure [8] shows the 
ratio of the magnetization to the exact value, computed 
using both methods. The relative difference reaches 10~^ 
close to the left boundary of the critical region window. 
We also checked that the low- and high-temperature sus- 
ceptibilities are not very sensitive to the way the mag- 
netization is calculated and Fig. [8] gives preference for 
using the definition p2p which we use in this paper for 
the calculation of the magnetization and the magnetic 
susceptibility. This was mentioned already in the pa- 
per [33|. 

Figure [S] shows a comparison of MC data for the mag- 
netization with the expansion (I21[) of the exact result, 
where the MC data for several lattice sizes are divided 
by the exact value. We see that the MC data coincide 
rather well with the exact result down to the reduced 
temperature \t\ sa 0.003. This value thus defines the left 
boundary of the critical region window in our subsequent 
analysis. 

The data for the effective amplitude Bs — ms/\T\^^^^ 
(see Fig. [T0| were fit with the expression 



Bs=Bo + Bi\t 



1 2/3 



Oi\r\) 



(45) 



which, up to this order, agrees with the expansion ()21|) . 
For L = 162, in the temperature window \t\ € [0.003 — 
0.012] (compare with Fig.|9]), the estimation gives an am- 
plitude of Bq = 1.2842(4) which is three standard devia- 
tions off the exact value Bq — 1.28326 . . .. For L — 243, 
estimated within the appropriate temperature window 
\t\ e [0.002 - 0.012], the value of the critical amplitude 
Bq — 1.2833(2) is in excellent agreement with the exact 
value. In both cases we found that in the temperature 
windows it is sufficient to just include the first correction- 
to-scaling term and we can neglect the second one in 



TABLE III: Results of fits to the MC data for the effective 
energy amplitude in tlie low-temperature phase. 



fit window 



£i(-) E2i-) 



E^i-) E4-) 



exact 



1.2783 . 



-0.8813... 0.5033. 



0.5401 . 



0.001<|t|<0.56 1.2771(10) -0.80(4) 
1.2946(4) -1.42(1) 



0.86(8) 0.99(5) 
-0.68(1) 



0.001<it|<0.32 1.2779(14) 
1.2834(8) 



-0.84(6) 0.86(8) 0.99(5) 
-1.13(2) -0.68(1) 



O.OOKltkO.l 



1.2795(13) 
1.2796(6) 



-0.99(4) 
-0.99(9) 



0.01(9) 




0.001<|t|<0.01 1.2795(17) -0.99(5) 
0.003<|t|<0.01 1.2779(23) -0.95(6) 



the fit. The estimated values of the first correction-to- 
scaling amplitude are i?i = —0.592(11) and —0.582(6) 
for L — 162 and L = 243, respectively, which are in 
good agreement with the exact value Bi = —0.58982 . . .. 
Figure [To] shows the fit to the effective amplitude of mag- 
netization data for L = 243. 

We also checked the equidistribution of magnetization 
moments (rnf) and polarization moments (p") with n — 
1, 2, 4 over the three sub-lattices Li and Ljk and found it 
valid within statistical accuracy. 



E. Polarization 

The critical amplitude Pq of the polarization 

can be estimated in the same manner as the magneti- 
zation critical amplitude Bq. The final value is Pq = 
1.2104(3), to be compared with the exact value [3J| 



Po = 2 



¥(ln(y2 + l))^V3^ 



1.20987. 



We may also estimate the ratio Bq/Pq ~ 1.061(1) which 
is close to the apparently exact value 3-\/2/4 « 1.06066. 



F. Low-temperature susceptibility 

Let us now come to the main subject of the present 
paper, the magnetic susceptibility amplitudes which are 
not known analytically. We can estimate them from our 
MC data using the same type of analysis we performed 
for the specific heat in the previous subsection. For an 
additional control of the analysis we compare our MC 
data with the available series expansions (SE) data. The 
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FIG. 8: Ratio of MC data with L = 363 for the magnetization 
to the exact result. The magnetization ttimc is computed 
using Eqs. pn|l (circles) and P^ (triangles). 



10 



available SE data [50| are short, however, and could not 
be used for a reliable estimation of amplitude values. 

Figure [TT] shows MC and SE data for the effective 
amplitude susceptibility F'^ = A:bTx-|t|^/^ in the low- 
temperature phase as function of — r. The solid line rep- 
resents a direct SE data summation, while the dashed 
line is the Fade approximant of SE data [51|. MC and 
SE data coincide well for \t\ > 0.01, and the discrepancy 
for smaller |r| is not surprising due to the short SE data. 

The parameters of the fit to the MC data shown in 
Fig. [TT] according to the expression 



effective amplitude data with the expression 



F' 



cff 



r- + ^r|r|^ 



D^\r\ 



D^\r\^ 



(46) 



are given in Table HVl Clearly, we can accept as the final 
and very conservative estimate the value F_ = 0.0681(1). 
Using the exact values for ^o (= ^o(— ) = ^o(+)) and 
Bq, we can estimate from this value the universal ratio 

i?c, = aAor-/B^ = 0.00517(7). (47) 

We also estimated the low-temperature critical ampli- 
fy) 
tude of the polarization susceptibility F_ in the same 

manner as for the magnetic susceptibility. The result is 

r^' = 0.061(1) and F_/F^^ « 1.11 which appears to be 

close to the ratio (Bo/Pq)^ = 1.12499 . . .. 



G. High-temperature susceptibility 

Figure [T^] shows MC and SE '5T\ data for the suscep- 
tibility x+ ill the high-temperature phase as function of 
T. The data coincide well at large enough reduced tem- 
perature T < 0.5 and diverge at small r because of the 
small number of terms in the SE available. 

The effective amplitude of the high-temperature sus- 
ceptibility is shown in Fig. [TS] Our results of fits of the 



rf ^ r^ 



D+T^ 



d;t 



D^T^ 



(48) 



are presented in Table [Vj 

The final estimate of the high-temperature suscepti- 
bility amplitude is F+ = 0.265(5). This implies for the 
universal susceptibility amplitude ratio the central esti- 
mate 



7^x^F+/F_= 3.9(1), 



(49) 



in very goo d agreement with the analytical predictions 
ofRefs. M and M. 



VI. DISCUSSION 

It is a widely accepted believe that there are four 
known models in the 4-state Potts model universality 
class. Besides the 4-state Potts model itself [lg|, these 
are the Baxter- Wu model [21[, the Ashkin- Teller model 
with some particular values of parameters [2J] , and the 
Debierre- Turban model J53{ with some particular value 
of parameter. Table IVII summarizes the known knowl- 
edge of universal amplitude ratios for the first three mod- 



c 



els, where we also included the universal ratio i?, 
aAo{+)T+/BQ. There are no estimations made for the 
Debierre- Turban model. Clearly, all estimates for F_|_/F_ 
from Monte Carlo simulations and series expansions are 
systematically smaller than the analytical predictions, 
also compatible within error bars, with a higher devi- 
ation reported in 1541 . There are some visible deviations 
of the result published in |20|[60|. The analysis of this 
quantity presented in Ref. [20| is based on the inclusion 
of logarithmic corrections to scaling (both multiplicative 
and additive) in the fit. This procedure is a bit risky 



1.002 
1.001 
1.000 
0.999 



.<" 0.998 
0.997 
0.996 
0.995 



i ? : I 



I * i ? ! i ^^ 



• L=162 
■ L=243 
^ L=363 



0.000 0.002 0.004 0.006 0.008 0.010 0.012 



1.275 

^ 1.270 

~ 1.265 

E 1.260 

1.255 

1.250 



0.002 0.004 0.006 0.008 0.010 0.012 



FIG. 9: Ratio of MC data for the magnetization nis to the 
exact result rrioxact for three values of lattice size. Error bars 
are only shown for the lattice size L = 243 (open squares). 



FIG. 10: Effective amplitude of the magnetization: Fit with 
expression (|45l) (sohd hue) to the data for L — 243 (solid 
circles) . 
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FIG. 11: Effective amplitude n of the magnetic susceptibil- 
ity in the low-temperature phase: MC data (circles), SE data 
(solid line), Fade approximant to SE data (dashed line). 
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FIG. 12: Magnetic susceptibility x+ in the high-temperature 
phase: MC data (circles), SE data (line). 



although one does not have to do something else. At 
the same time, it is argued in [20| that the universal 
combination R^ should not contain any logarithmic cor- 
rections in the effective estimation through the function 
Rc{\t\)- = a{a - l)(e-(|r|) - eo)x-/m2(|r|) and, in- 
deed, there is a good coincidence of results reported in 
Refs. [20, [2I] , and in the present paper for i?^ . By ana- 
lyzing the data in Table IVTl we may conclude that there 
are definite overestimations of the critical amplitude T ^ 
in Ref . [20| . This is possibly due to the large background 
terms (nonsingular contribution) in the high-temperature 
susceptibility. We have to note that there is only one di- 
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TABLE IV: Results of fits to the MC data for the effective 
low-temperature susceptibility amplitude Fl . 



fit window 


F_ 


Dr 


D2 


D^ 


0.002<|t|<0.5 


0.06814(5) 


-0.423(3) 


0.491(12) 


-0.11(1) 


0.002<|t|<Q.25 


0.06751(6) 
0.06800(3) 


-0.355(6) 
-0.406(2) 


0.170(24) 
0.401(3) 


0.25(3) 



0.003<|r|<0.25 


0.06931(9) 
0.06893(5) 


-0.468(7) 
-0.435(2) 


0.592(29) 
0.449(3) 


-0.15(3) 



0.002<|t|<0.05 


0.06804(31) 
0.06909(10) 


-0.308(38) 
-0.444(7) 


-0.32(22) 
0.473(20) 


1.26(34) 



0.002<[r[<0.025 


0.06655(8) 
0.06669(3) 


-0.273(7) 
-0.286(1) 


-0.042(23) 







0.003<|t|<0.025 


0.06879(13) 
0.06715(4) 


-0.420(9) 
-0.295(1) 


0.39(3) 







0.003<|t|<0.011 


0.06770(6) 


-0.313(2) 








0.004<|t|<0.011 


0.06801(9) 


-0.32(3) 









FIG. 13: Effective amplitude of the magnetic susceptibility in 
the high-temperature phase. 



rect estimate of the universal ratio r7^/r_ published in 
[20|, which is not consistent with analytical predictions. 
More work should be done to clarify this issue. 

Finally, we may conclude that our analysis of criti- 
cal amplitudes of the Baxter- Wu model produces uni- 
versal amplitude ratios consistent with the analytical re- 
sults obtained by Cardy and Delfino for the 4-state Potts 
model 28] and by Delfino and Grinza for the special case 
of the Ashkin-Tellcr model [30]. 
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TABLE V: Results of fits to the MC data for the effective 
high-temperature susceptibihty amplitude F^. 



fit window 


r+ 


D+ D+ D+ 


0.002<|r|<0.77 


0.265(2) 
0.276(1) 


0.79(7) -0.93(18) 0.79(12) 
0.34(2) 0.25(2) 


0.002<|r|<0.43 


0.267(3) 
0.272(1) 


0.68(12) -0.59(34) 0.53(27) 
0.45(3) 0.09(4) 


0.002<|r|<0.025 


0.266(3) 
0.265(1) 


0.74(48) -0.7(1.5) 
0.549(3) 


0.004<|r|<0.025 


0.270(3) 


0.52(6) 


0.002<|r|<0.01 


0.268(4) 


0.57(12) 


0.004<|r|<0.01 


0.269(7) 


0.54(18) 



TABLE VI: Universal combinations of critical amplitudes for 
the two-dimensional models in the 4-state Potts model uni- 
versality class. 



model 


approach 


A+/A. r+/r_ Ft/F- Rc Ri 


Ref. 


4-state Potts model 


analytical 

MC 

SE 

MC and SE 


1.0 4.013 0.129 0.00508 0.0204 
3.14(70) - 0.0068(9) 0.021(5) 
3.5(4) _ _ - 

1.000(5) 6.49(44) 0.154(12) 0.0052(2) 0.0338(9) 


[27,^ 
[54] 
[55j 
[20] 


Ashkin- Teller model 


analytical 


4.02 0.129 


[30] 


Baxter- Wu model 


MC 


0.995(5) 3.9(1) - 0.00517(7) 0.0201(5) 


present 
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We have to note some recent discussion of the weak uni- 
versality of correlation length amplitudes [lOl. [Till and as- 
sociated quantities like Binder cumulants [la . [la | in the 
case of anisotropic models. 

At least in the directions of the "space" dimensions, as is 
the case for the transfer-matrix method where the "time" 
direction can be taken to infinity. 

The temperature at which thermodynamic quantities, 
such as the specific heat and susceptibility, reach a max- 
imum value. The thus defined pseudo-critical tempera- 
tures reach Tc in the limit of infinite lattice size, and the 
maximum value of thermodynamic quantities diverges 
accordingly. 

It has been demonstrated that such implementation is 
very helpful in simulations of the bond-percolation prob- 
lem on the square lattice [441 ]. 

We mistakenly mentioned in [20] a wrong estimation of 
the value of F+ZF- for the Baxter- Wu model. 



